NPS ARCHIVE 
1964 

HILDER, L. 



EXTENSION OF THE ROOT LOCUS 
TECHNIQUE INTO THREE DIMENSIONS 

LEONARD 0. HlLDER, JR. 

AND 

RICHARD G. GEER 




Library 

XJ. S. Naval Postgraduate School 
Monterey, California 



EXTENSION OF THE ROOT-LOCUS 



TECHNIQUE INTO THREE DIMENSIONS 



by 

Leonard 0. Hi.lder, Jr. 

'/ 

Lieutenant Commander, United States Navy 

and 

Richard G. Geer 

Lieutenant, United States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 



MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 



19 6 4 



feLTeJ^ 



Lfh’-n-r 

U. S. N mil Po tg'fuluate Schooj 
Monte; ey, California 

EXTENSION OF THE ROOT -LOCUS 
TECHNIQUE INTO THREE DIMENSIONS 
by 

Leonard 0. Hilder, Jr. 
and 

Richard G. Geer 

This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 
from the 

United States Naval Postgraduate School 



ABSTRACT 



The root- locus technique used in the determination of system response 
has been a powerful tool for the engineer. Shortcomings of the technique, 
however, particularly in regard to the locations of imaginary roots and 
the frequency stability of the system for changes in the system gain, 
have made a supplementary technique desirable. Such a technique utiliz- 
ing the projection onto the K- plane of the three-dimensional root- 

locus has been developed. The use of the root-locus plot in conjunction 
with the K- ^<-0 trace eliminates these shortcomings and provides the 
engineer with still another tool in his quest for a more complete method 
of systems analysis. 
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INTRODUCTION 

For many years engineers have been using rooe-loeus methods of 
solving high-order polynominal characteristic equations of control 
systems. As an example, consider the simple servo system shown in 
Fig. 1.1. 




y ( a ) 



In general, the blocks G and H represent any pole-zero combinations 
which might describe the system. The characteristic equation of the 
system is 

\ + GH = O , 



which can be written in general terms as 

• (4+tv,) (a +?„*) .. i i 

where n • 1, 2, 3, - - - 

K • root- locus gain 

P * pole of G 
G 

P “ pole of H 

H i 

Z_ * zero of G 

W 

Z » zero of H 
H 

A • 

s ■ Laplacian variable = CT~ 4 
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Obviously, equation 1.1 could be a polynorainal of any order, and the 
analytic solution of this equation becomes tedious and impractical as 
n Increases. Conaaquantly , the root-locua method of solution based upon 
the pole-zero configuration of the open-loop transfer function GH, which 

i 

is both simple and straightforward, has been widely used in engineering 
practice. Such a plot provides a wealth of information about the system 
in question and is relatively easy and quick to obtain. There are, how- 
ever, certain shortcomings to the method. Notably, the root-locus plot 
of a given system response does not permit ready evaluation of the loca- 
tions of roots of the characteristic equation for a given value of the 
, root-locus gain, K, or over a specific range of K. A trial and error 
measurement procedure is required, and for very small K where the roots 
lie close to the poles, this graphical method usually fails completely. 

In addition, there is no information concerning sensitivity of the system 
to changes in gain. 

The first of these shortcomings has been partially resolved by Dr. 

R. C. H. Wheeler, of the U. S. Naval Postgraduate School. Substituting 
A = <T44 v> into equation 1.1 and letting ^U)=© , Dr. Wheeler was able to 
obtain a plot of K vs. f (called the "Wheeler plot"). With the Wheeler 
plot, the real roots of the characteristic equation can be determined, 
if they exist, for any value of the root-locus gain. The most valuable 
contribution of the Wheeler plot was the prediction of a maximum or a 
minimum value of K for those at which .the root-locus emerges from or 
enters the (F-axis. thus, the value of K at these points, called 
and respectively, may be quickly determined from the root-locus 

* j 

plot without the necessity for plotting the Wheeler plot. An example 
of the Wheeler plot is shown in Fig. 1.2 for. the open-loop transfer 
function, GH * K U+l) . 

4-U + O 
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Root-uocvjs 




ROOT-LOCUS AND WHEELER PLOTS FOR GH - FOR 

4 . ( 4 + 1 ) 

ARROWS DENOTE INCREASING K. 

Fig. 1.2 



K 2.0. 
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However, no information is available frcxn the Wheeler plot concerning 
complex roots, nor even about the real parts of such complex roots. 

Returning to equation 1.1 as the first step in searching for 
answers to these shortcomings, the substitution of 4- <T-»gu> immediately 
revealed this to be an equation in the three variables, (T , , and 

K. Since the root-locus is a plot of 0" vs. for varying K, such a 

plot is actually the projection upon the plane of the three-dimen - 

sional plot of equation 1.1.. Such a three-dimensional plot for the pre- 
vious example, GH ■ “4“^ » * 8 shown in Fig. 1.3 and is interpreted 

4 1 A+i; 

as follows: 

1) Curves B and C represent the three-dimensional plot of 
the characteristic equation. 

2) Curve B, for which is the Wheeler plot for the 

system. This curve is identical to the plot shown in 
Fig. 1.2 for this same system. 

3) Curve A, being the projection of the three-dimensional 

curve upon the plane, is the root-locus plot for 

the system. 

4) Curve D, the projection upon the plane, is called 

the trace 1 ' in the succeeding development. 

Curve D of Fig. 1.3 showed great promise of providing a solution to 
the shortcomings mentioned previously connected with the root-locus plot: 
the imaginary part of any complex root is immediately available for any 
value of K, and the frequency stability for varying K is obvious by in- 
spection over any given range of K. The problem, then, became that of 
determining whether the K trace for any given system could be 
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t 



44 *° 



K 




THREE-DIMENSIONAL PLOT OF 1 + GH - 0 FOR GH 



k <■»-+«) 

A(*+l) 

FOR K^O AND FOR >0*0 . ARROWS DENOTE 



INCREASING R. 
Fig. 1.3 
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e 

• ( i, 

v 

sketched quickly and accurately. Speed and ease of sketching were 

i 

« 

paramount requirements; otherwise, a graphical trial and error pro- 
cedure utilising the root-locus would provide the required Information 
without the necessity for making an additional plot. 

Two methods of attack were pursued in seeking a solution to the 
problem of rapidly and accurately sketching the K — 4 VO trace for any 
given open-loop transfer function: 

METHOD I : The location of poles and zeros on the K-axis 

through some suitable transformation upon the poles and 

I 

zeros of the given open-loop transfer function. 

In searching for a suitable transformation to use, an independent 
variable, M, in the K plane was defined as follows: M A K-M<q 

where M was a vector in the K - 4 U 0 plane from the origin to any point 
(K, ) on the K -,ju> trace. M, then, corresponded to the Laplacian 

vector, s, in the <P> 4 U> plane. These vectors are depicted in Fig. 1.4. 




Fig. 1.4 
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The attempt was then made to express the open-loop transfer function 
in the form 



GH a- - . _ 

where small letters denote zeros on the K-axis, capital letters 
(except M) denote poles on the K-axis, and n represents the system 
"type" when the transformation has been completed. This method broke 
down immediately, however, as soon as it was attempted with the simplest 
second-order system, GH ” 4r • Substituting into the equation 

1 + GH • 0 resulted in •> 

+ ( 4<*>) X + K 4 ^ 2 (TW = O 1.3 

Inspection of equation 1.3 revealed the impracticability of separating 

M and obtaining the desired form of GH. This can be seen more clearly 

o 

by dividing equation 1.3 by K + (^U)) , resulting in 

G <ri£±i2^ = _, 

Equation 1.4, then, shows the two difficulties encountered using this 
method: 1) the inability to form the vector M, and 2) the inability 

to completely separate out the independent variable <T . This method 
was' attempted on several different transfer functions with the same 
(and often worse) result. Consequently, method I for obtaining the K-^ w 

trace was abandoned. 

( 

METHOD II : The location of critical points on the K-axis such 

that a set of rules similar to those used for making a root-locus 
plot might be established for sketching the K trace. 
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Basically, this method involved looking at the projection of the 
three-dimensional curve of the characteristic equation upon the 
plane. Prom this projection, then, critical points were established 
as appeared most promising for the system under consideration. Al- 
though the results in some cases were slightly different than antici- 
pated, the use of this method provided a successful solution to the 
problem of obtaining the K- 4 ^ trace quickly and accurately. For 
the case of characteristic equations in the form of second-order poly- 
nomials, the trace was determined ^analytically - the develop- 

ment for this case is traced through in chapter II. For characteristic 
equations in the form of third or higher order polynomials, an approxi- 
m ate graphical method of obtaining the K -4UJ trace was developed as 
shown in chapter III. 

Utilizing method II, therefore, it was shown that the K -4W trace 
for any given open-loop transfer function can be sketched quickly and 
accurately using the concepts developed in the succeeding two chapters. 
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SECOND-ORDER SYSTEMS 

2.1 DISCUSSION . For the purpose of this development, a second-order 
system is defined as a system having a characteristic equation in the 
form of a polynomial of degree less than or equal to two. For such a 
system, a relatively easy analytic solution to both the root-locus plot 
and the K—Jco trace was obtained. 

The Laplacian variable,- s - , is a vector whose variable com- 

ponents (in the mathematical sense) are <T and 60 , with ^ = *$-1 being a 
90° rotation of the <*> -axis from the T -axis. In engineering practice, 
however, the variables are usually taken to be 9* and jcO — hence, 
the term "K- jw trace". Although this term is used throughout this 
paper, the second variable will be <0 rather than J(0 for purposes 
of the mathematical development in this chapter, although the variable 
plotted will be joJ . Elsewhere in this paper, j«*> will be used as the 
variable in accordance with accepted engineering practice. The illustra- 
tions shown in this chapter will be plotted only for positive values of 
K, since negative values of K result in a 0° root-locus. The K-JfW 
trace for negative values of K is discussed in chapter III. 

As a general second-order system, the open-loop transfer function, 

GH » was used for this development. 

2.2 ROOT-LOCUS PLOt I The characteristic equation, 1 + GH ® 0, of the 
general second-order system became; 

(b+A*W * F+OCsO. 

> 

1 

Although this, 'is not a new concept, the work was original and 
necessary for the development of the K- jto trace following in 
section 2.3. 
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Substituting s ■ *T+Ju> resulted in the equation, 

-0>+Aic)(i) J +(fe+BKV+ F-KLk] +j<o[2(t>+AKV +E+Bk] = O. 

Noting that this equation was in the vector form, » it 

followed that 0 4 - ■ 0 and ■ 0. Writing these expressions explicitly 
yielded the pair of equations, 

(6+akV 2 -(fe4A*)c»* + (E46K)r-v F4CK = O 2 .i‘ 

and 

doUfb4AK)r+ E+Bk] * O, 2 2 

The reason for using <0 as the third variable now becomes evident: 
equations 2.1 and 2.2 are the parametric form of the characteristic 
equation and therefore are the parametric form of the curve of the 
characteristic equation. These equations, of course, were obtained only 
under the assumption that the characteristic equation could be written 
in vector form, which necessitated using oJ as the variable. Equation 
2.1, with U> ■ 0, is the equation of the Wheeler plot. For this condi- 
tion (i.e., - 0) the root-locus is easily determined as well as the 

Wheeler plot, and no further investigation of this condition was deemed 
necessary. 

For the condition <*)/ 0 , however, the root-locus is not so easily 
determined, and the Wheeler plot is not defined. Now, the root-locus plot 
becomes a projection of the three-dimensional curve of the characteristic 
equation upon the plane. From equation 2.2 it was evident that UilO 

implied the variables and K could be determined in terms of each other, 

yielding: 
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Ibr + E 
£.At +• 6 



2.3 



and 



ij— _ i E + 6< 

*(D+AK) 



2.4 



The projection of the three-dimensional curve upon the T"-JcO plane, 
then, was obtained by substituting equation 2.3 into equation 2.1. 
Performing this substitution, collecting terms and completing the square 



in q* resulted in the equation, 



L AF-Ct>7 


* L CE-8F | 


f AF-Cbl* 




L T + AE-&& J 


+ " * aT^d + | 


a 

CD 

1 

tu 

/ 


2.5 



Clearly, equation 2.5 is the equation of a circle centered at 

_ AF -CD 

T a ■ AE-6D J <*> = O 

with radius 

It was concluded that the portion of the root-locus plot for which 
(O V O for all second-order systems is a circle. This circle can be 
determined exactly by using equation 2.5 in conjunction with the coef- 
ficients of the open-loop transfer function with two exceptions: 

1) For an open-loop transfer function containing no zeros, 

A ■ B * 0, and equation 2.5 no longer applies.' Returning 
to equations 2.1 and 2.2 in this case, the circle is de- 
generate into a straight line given by 

2) For an open-loop transfer function having a pole-zero 
configuration resulting in poles separated by zero (s) , 
or vice-versa, there is no root-locus for O . 
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Open-loop transfer functions with pole-zero configurations resulting in 
the applicability of the foregoing two exceptions present no difficulty 
in quickly sketching an accurate root-locus plot. The root-locus plot 

for all other second-order systems may also be quickly and accurately 

1 1 

sketched by substituting the coefficients of the numerator and denomi- 
nator polynomials of the open-loop transfer function into equation 2.3 
to determine the center and radius of the required circle. 

2.3 K- jtJ TRACE . * The K- jw trace was developed in two parts: 

the portion for which W « 0, and the remainder for which < 0^0 

For io=0, the K-jco trace must lie entirely on the K-axis. 

Since the Wheeler plot is the K- VT plot for lO = 0, the K - J OO trace 
for (0=0 must be simply the projection of the Wheeler plot into the K- 
axis. From a knowledge of the Wheeler plot, then, the K- Jo) trace for 
co m o can be determined from the real pole-zero configuration of the 
open- loop transfer function as shbwn in Table 2.1. 

Following the same procedure for 6 O which was used for the root- 
locus development in section 2.2, equation 2.4 Was substituted into the 
equation e>< ■ 0, resulting In: 

' ' 2.6 

*When the term "K" is used in conjunction with the "K- Jo) trace", 
it is inferred to be the root-locus gain of the open-loop transfer func- 
tion. .In general, this will also be the gain of = because 

A-D.-l, 

if both coefficients are present in any given open-loop transfer function. 
These coefficients are utilized in this development solely to allow the 
second-order terms in the open-loop transfer function to go to zero in 
special cases. 
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REAL POLE -ZERO CONFIGURATION 



K- \( j 3 TRACE ON K-AXIS 



I Single pole/zero. 

— . -A 






Entire K-axis. 



II Pole-zero combination (with 
root-locus on the T' -axis between 
the pole and the zero). 



AJU) 






^ K 



4k- 



-* <r 



HI. Pole-Pole combination (with root 
locus on the T -axis between the poles) 



Kemerc 



-*= 






K-axis from the origin to 
K EMERG 



4«VA 



-><r 



-> K. 






IV. Zero-Zero combination (with root- K-axis from K^^ to infinity, 
locus on the T -axis between the 
zeros) 



Kturr 









-> T 



^•K 



vevr 



V Double pole/zero. 



U0_ 



A 



None. 



r 



REAL PART OF K- Jco TRACE FOR CERTAIN OPEN-LOOP POLE-ZERO 

CONFIGURATIONS 



Table 2.1 
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Immediately, difficulty in obtaining a simple, analytic solution was 



encountered, because equation 2.6 is NOT the equation of a simple conic 



suggests a quartic equation. An analytic solution of the quartic equa- 
tion 2.6, although possible, would not yeild an "easy" method of deter- 

l 

mining the K-J&> trace. 



considered: (1) the "partitioning" method, and (2) an "approximation" 

method. 

(1) " Partitioning" method . The method has been generally 
described elsewhere as follows: any polynomial (up to 

about fifth degree) can be "partitioned" by dividing the 
entire polynomial by one or more terms of the polynomial, 

* yielding an equation in root-locus form, i.e., GH »-l, where 

GH is some combination of poles and zeros multiplied by some 
variable. For example, the polynomial 



section. 



•> 2 o 2 

The term 8ADKW 4 suggests a cubic, and the term 4A K 4 co 



Two possible graphical methods of solution of equation 2.6 were 



+ e a O 



4 I 3 

can be partitioned by dividing by ( % + b*. ), yielding 



the root-locus 




% 3 (■*-*-+ b') 



where c is the "root-locus" variable. In order to use 



this method, the original characteristic equation was 



utilized, and (Jw ).was considered as the third variable. 



yielding 



4^+BKv+F+CK+^lofj«j)4 2AKT(j<o) +E(jw) 
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Difficulty was immediately encountered due to the 
presence of the cross-variable terms, and no amount 
of juggling would enable the "root-locus 11 variable, <T , 



to separated in* such a manner that the "root-locus 
form" of equation 2.7 could be written. Basically, this 
method required the defining of a new vector , 

to provide a correspondence to the Laplacian variable in 
the root-locus plot, as was sh^>wn in chapter I. In- 
spection of equation 2.7 shows why this method could not 
be adopted - the required separation of variables could 
not be achieved. 

(2) " Approximation" method . Returning to the general 
second-order open-loop transfer function, 

6H " B-i+F 

the denominator waa divided into the numerator, resulting 
in 



KA + Vc(B^+C') 



<SH= -5- 






2.8 



where and c/= C--^ • Defining 

. S > L fc(B!Uc') j 

(<SH) * , it was noted that (GH) had the 

appearance of an open-loop transfer function similar to 

GH, the only difference being the number and positions of the 

zeros. Furthermore, 1 + (GH)* • 0 was known to yield the 

equation of an ellipse in the plane? This was proved 

by returning to equation 2.6; 

The two cross-variable terms both involved 

the coefficient A. -Since the presence of 
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these terms precluded equation 2.6 from repre- 
senting a conic section, their elimination (by 
setting A ■ 0) would result in the desired conic 
section. Performing this operation, collecting 
terns and completing the square in K resulted in 
the equation 

r <ci>-»e ~ | 3 v 

l K + <** 

^ (B*F-BcE +c i D) £ (B*F-Bce +c a 



2.9 



By inspection, equation 2.9 is the equation of 
an ellipse centered at [k - — = O J 

with semi -axes given by .1 and 

^ t> 

nJ t> (& l F - Bee +e 2 t>) , 

Since (GH)* yielded an ellipse in the K-^OJ plane, the 
relationship between GH and (GH)* was determined in an 
effort to obtain .information about the K-Jio trace result- 
ing from GH. This was accomplished by returning to the 
root-locus plots for the two transfer functions, since 
both plots must yield circles: 

The root-locus plot of GH was a circle having a 
center and radius as defined in equation 2.5. Plugging 
these formulae with A ■ 0, B ■ B / and C ■ C / resulted 



in the center and radius of the circle for (GH) as 

i 

c 

[T 6 * b' * t * >s ° 1 



Substituting B 



_ AE 

B ~ "BT 






and C 



' - C--^ , 



it was 



found that the center and radius of the circle for 



(GH) expressed in terms of the variables of GH were 
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identical to the center and radius of the circle for 
GH! This is shown graphically in Fig. 2.1. The basic 
relationship between GH and (Gh/ can be seen by in- 
spection of Fig. 2.1. The effect of the division per- 
formed to obtain (Gh/ was to move the zero at to 

the center of the circle and the zero at to in- 

finity. The root-loci of the two functions for ua^O 
are identical in shape and position, the only difference 
being the values of the root-locus gains at correspond- 
ing points on the circles. 

With this correspondence established, it was inferred 

that the K- Jw trace for GH must be in the shape of some 

sort of "stretched" ellipse with the same maximum L) as 

that for the true ellipse of (GH) 7 . To quickly sketch 

the K-jCO trace for GH, then, it would be sufficient to 

A 

determine the key points given by = the value of 

the root-locus gain for which the root-locus departs the 

A 

T" -axis, Kgty<j< — the value of the root-locus gain for 
which the root-locus enters the T -axis, and — the 

value of the root-locus gain at the peak (maximum 63 ) of 
the root-locus circle. Since and Kg^ must both 

lie on the K-axis, and since the value of CO corresponding 
to is known from the radius of thq root-locus circle, 

the key points of the K-jiO trace are determined, 
and Kgtyj are quickly available from the Wheeler plot, and 

is equally quickly available by a simple measurement 
on the root-locus plot. With these key points, then, and 
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ROOT-LOCUS PLOTS FOR GH AND (Gh/ FOR GH 
AND 

Fig. 2.1 



& W~ ) k>/ ° 



} 



1 

I 
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/ 



knowing the general shape of the curve, the K-J<o trace 
can )>e quickly sketched. 



For the special case in which OH contains no zeros, neither aquation 
2.9 nor the approximation method provide a valid solution to the K - Jd 
’ trace. This was intuitively obvious, because an open-loop transfer 

M 

function containing no zeros yields a root-locus plot for d/O which 
is an OPEN curve, i.e., the curve does not close on itself except at 
infinity. Consequently, the projection of the three-dimensional curve 
onto the K-^’w plane should also be an open curve. Further, it must be 
a conic section as may be seen by setting A ■ B « 0 in equation 2.6. 

The K-jco trace for this case, then, must be a- straight line, a para- 
bola or a hyperbola. Substituting A « B » 0 into equation 2.6 and re- 
arranging terms resulted ip the equation 



Equation 2.10 is the equation of a parabola symmetrical with respect to 



l 




2.10 



the K-axis with vertex at l ~ 3 ' » 40 “ 0], thus substantiating 

the intuitive result anticipated. 



It was concluded that the K- traces for all second-order systems 
may be obtained quickly and accurately utilizing Table 2.1 for the 
portions for which O- 0 and Table 2.2 for the portions for which **>^0. 
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OPEN -LOOP 
TRANSFER FUNCTION 


SHAPE OF 
K-jw TRACE 


METHOD OF 
OBTAINING 


ALTERNATING POLES 
AND ZEROS 


STRAIGHT LINE 
(K-AXIS) 


PROJECTION OF 
WHEELER PLOT 


NO ZEROS 


PARABOLA 


EQUATION 2.10 


1 ZERO 


1 ELLIPSE 


EQUATION 2.9 


2 ZEROS 


"STRETCHED" 

ELLIPSE 


APPROXIMATION 

METHOD 



SHAPE OF K-Joj TRACE FOR VARIOUS OPEN-LOOP POLE-ZERO 
CONFIGURATIONS, oa^O. 

Table 2.2 



2.4 EXAMPLES . Four examples were selected to illustrate the various 
solutions for the root-locus plot and the K-jto trace for second-order 
systems. 



1 . 



<3H = 






This is the case of alternating poles and zeros, and the plots are 
shown in Fig. 2.2. 




The root-locus plot for (&i-0 is a straight line parallel to the 
^09 -axis and passing through the point The K- trace 

is a parabola symmetrical with respect to the K-axis with vertex 
at [ K m £-Z m “ 0]. These curves are shown in Fig. 2.3. 

3 . 6H s B=b- 1 F-Z 

4*4344 A ** U '> ^ F Z ' 



The root-locus plot is a cir cle centered at ( V" • -5, & ■ 0) 
with radius + a * Z.AiA, 
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10-3 

The K- trace 18 an ellipse centered at ( K = — ; — = 7, 

w • 0) with (semi-axis)^, - IS +2S m ■ 6.928 and 

(aemi-axis)^^ ■ ^(^F ■ 3.464. These are plotted In Fig. 2.4. 

. ... K{*+3)f4+4) K(4*+l4+U) . . ^ , 

4 - = ?»+fffi+ay = P=i. 



l _-\2 

The root-locus plot is a circle centered at ( V *= — 3 ~~] ~ -2.5, 

(*> « 0) with radius R " IW + (i.£> 1 - ^^a7? » o. U7. 

The K-joo trace is a "stretched" ellipse with key points Kgj.jgj^g 

0.0718, K ■ 13.92, Kuav " 1*0* These curves are plotted in 
ENT 

Fig. 2.5. 
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ROOT-LOCUS PLOT AND K- jUJ TRACE, 
GH - *JL*tQ.l**$ 

A.(*+2) 




<■ 

3 
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Fig: 2.3 



iw trace, 
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ROOT-LOCUS PLOT AND K - \U> TRACE, 




!<*>♦ 



t 
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HIGHER-ORDER SYSTEMS 



3.1. PROCEDURE . The methods developed in the previous chapter could 
not be readily extended to give solutions for higher order systems. 
Therefore a different approach to the problem was indicated. 



Since the desired procedure was to be an addition to the root- 
locus method it seemed logical to try to extend root-locus techniques 
to achieve the desired projection. On the s-plane, the root-locus 
must satisfy both the magnitude and angle criterions. Those points 
which satisfy the angle requirement can be accurately determined by 



guide lines were not available. In practice, the centroid is computed 
and the asymptotes are drawn, then the root-locus can be rapidly "rough- 
ed in" with the Spi-rule used mainly as a check and to provide added 
accuracy if desired. Thus, if a method of determining the centroid 
and the angles of asymptotes could be found in the plane, per- 

haps with practice the desired K- ,jtt> projection could also be "roughed 
in" and then checked for accuracy as necessary by some other means. 

3.2 ASYMPTOTES . To develop a method of finding asymptotes in the K-ju> 
plane it was helpful to consider first why there were asymptotes in 
the <T-jioplane. As the parameter, K, is allowed to become very large, 
those roots which are approaching zeros at infinity are distant from 
the open loop poles and zeros. The angle condition must be met by the 
contribution from the excess poles since the algebraic sum of the angles 
from the remaining poles and zeros will effectively cancel. Thus, as K 
goes to Infinity, the root locations allowable are limited to lying on 




would be a very tedious procedure if other 
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asymptotes having angles which, when multiplied by the number of excess 
poles, satisfy the angle criterion. The angle of the asymptote must be 
measured at the effective center or centroid of the open loop poles. 



In much the same manner the value of gain required for a point on 
the s-plane asymptote which is distant from the open loop poles and 
zeros is determined by the product of the distances from that point to 
the excess poles, since the distances from the other poles and zeros will 
effectively cancel. An approximation for the gain at this point is the 
distance to the system centroid raised to the power n, where n is the 
number of excess poles. Referring to Fig. 3.1, the distance from a 

point, R, on the s-plane asymptote to the centroid is 

J _ to, 



3 SW> 

where U), is the value of u> at k and o L is the s-plane asymptote angle. 
Thus the magnitude of the gain is given approximately by 



the approximation being valid if U), is very large. Consequently, 
solving equation 3.1 for U>, and assuming U), very large gave 

CO a (S ivv •O’VTkT 3.2 



where ■ the angle of the s-plane asymptote, and n * the number of 
excess poles. Equation 3.2 was obviously not linear in K and UJ , 
implying the asymptotes in the s-plane did not transfer into straight 
line asymptotes in the K- jU) plane under the approximation used. 

This implication appeared to be a serious drawback to this method, 
because curved asymptotes could not, in general, be quickly and easily 
plotted. However, a compromise solution was reached: by letting the 
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not* i 



8-plane Configuration 




s.^ * 



gain variable be wr , equation 3,2 Is the equation pf a straight line, 
since sin is a constant* for any given asymptote. Further, this 
straight line, which is the s-plane asymptote transferred into the ^VlKl 
plane , always passes through the origin with slope ■ sin«c, . Consequent- 
ly, the only knowledge of the system required to quickly plot the asympto- 
tes is n, the number of excess poles of the system. The advantage of 
having straight-line asymptotes, then, made the adoption of the ^f|K| -jw 
plane advisable for obtaining the required trace. The term "K- 4 <*> trace" 
will still be utilized in this discussion, but the abscissa is implied 
to be V/IM * A discussion of the relative advantages and/or dis- 
advantages of adopting this procedure is contained in section 3.5. 

Since the root locus is symmetric about the real axis the K-ju> 
trace will be symmetric about the axis. Thus, as in the root lotus, 

only the upper or positive ^t*>- axis need be used since the asymptotes 



and trace for the negative half would be the mirror image of the posi- 
tive half. 

It should be noted that for negative values of K the root locus 
obeys the same magnitude criterion as for positive K. Therefore the 

f 

asymptote in the W plane will obey the same equation for both 

positive and negative K, that is, equation 3.2. Thus the trace for 
negative K could be plotted on the same axes as for the positive K trace, 
but for clarity the following convention was established: The positive 

trace is drawn on the + *J|K| coordinate and negative trace on the 



It should be noted that for higher order systems there will be n 
asymptotes on the s-plane, each of which will be defined by an angle oC . 
There will be a K-4<*> asymptote for each s-plane asymptote defined 
by equation 3.2 utilizing the appropriate tc • 
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— ’ty iKl axis, with the 4 W -axis passing through the zero point of the 
abscissa. Care must be used in plotting the K- ju> asymptotes, because 
they are not , in general, equally spaced (in the angi^lar sense) as are 
the s-plane asymptotes. T^e K-j 10 asymptotes for positive K and those 
for negative K must be sketched separately, for the angles of the s-plane 
asymptotes are different for positive and negative values of K. This is 
shown graphically in Fig. 3.2. 




Negative K 



ROOT LOCUS AND K-4 W ASYMPTOTES FOR GH « FOR POSITIVE AND 

NEGATIVE VALDES OF K , ^tO 2 : 0. 

Fig. 3.2 
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3.3 K-sju> TRACE . Having determined the asymptotes in the 
plane, the first step in determining the K- 4 U > trace involved a review 
of the procedures for sketching the K> jtu trace of second-order systems 
as outlined in Tables 2.1 and 2.2. 

As a first consideration, the real double pole/zero combination 
(as shown in category V of Table 2.1) produced no trace on the K-axis, 
and therefore had no trace on the wr -axis, since the latter axis 
differed from the former only by a scale factor. Using Table 2.2 (for 
no zeros) and equation 2.10 yielded the equation 

Ol ) 1 = K 3.3 

Transforming to the 4iFJ coordinate, equation 3.3 became 

u>= Vk 

which is precisely the equation of the asymptote in the YlKl-jU) pi ane. 

I < 

Therefore, the K- jto trace for such a system lay entirely on the asymp- 
totes, as shown in Fig. 3.3. 




Fig. 3.3 
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The extension to a real multiple pole of any order, then, was simply 
a combination of the foregoing result and successive applications of 
Table 2.1, resulting in the conclusion that a real multiple pole/multl- 
ple zero of any multiplicity will have a K-juj trace which coincides 
precisely with the K- ^co asymptotes. For poles of odd multiplicity, 

rl , 

of course, the entire ^/| K| -axis will be part of the trace, since 
the odd pole falls under category I in Table 2.1. 

As the poles of a real .double pole combination are separated along 
the o- -axis, the change in the R- .jU> trace (from that obtained for 
the double pole) was expected to be most noticeable near the origin, 
since the curve must become asymptotic for increasing values of K. 
Further, the value of K at which the trace departs the K-axis has al- 
ready been determined to be (which becomes ^Kem£ R (j on VF 

-axis). Consequently, the K-ju> trace was sketched as shown in Fig. 
3.4. 




ROOT LOCUS PLOT AND K - ^ 



Fig. 3.4. 



,K*C> AND 
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For this special case, i.e., an open-loop transfer function GH “ (4 + c «j(l>+b) * 

an interesting property of the K- 4 W trace became apparent. The K-jto 

trace in Fig. 3.4 bore a striking resemblence to the root locus plot of 
k 



GH 



(44\aV 



v , as may be seen by comparing Fig. 3.4 with Fig. 3.5, 



if the pole, a, was relate^ to NK 



EMERG 




Since all points on the root locus shown in Fig. 3.5 must necessarily 
satisfy the angle criterion, it was felt that a similar angle criterion 
might be established for the K-jtu trace. Locating a simulated double 
pole at and another simulated double pole at ~ w^^lERG in Fig * 

3.4, the K- \U> trace satisfied a 180° angle criterion with respect to 
these simulated poles - precisely the same angle criterion established 
for the root locus plot; It was concluded, therefore, that the K-4 60 
trace for an open-loop transfer function consisting of two real poles 
could be sketched quickly and could be checked directly with a Spi-rule 
(utilizing the simulated poles described earlier) if greater accuracy 



was; required. 
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For an open-loop transfer function consisting of a pair of com- 
plex poles, the K- trace was easier to sketch than for any other 

polt-caro configuration. Sinca the jU) -axis was common to both tha 
root locus plot and the K-njU> trace, and since the trace must be 

• M 

asymptotic, it was a simple matter to sketch the trace departing the 

h 



j Co -axis at - the complex'''values of the poles and approaching the 45° 
asymptote. The K-JCO trace for such a pair of complex poles is shown 




■> 4TkT 



ROOT LOCUS PLOT AND K- 4<*> TRACE FOR GH - f- , 

Kle AND AWJro 
Fig. 3.6 

The procedures to use in sketching the K- A»*> traces for the fore- 
going "basic pole configurations" of an open-loop transfer function are 
summarized in Table 3.1. These procedures are equally applicable by re- 
placing the word "pole" by the word "zero", " Aj K^^," by. " and 

"K - 0" by "K « infinity" where they appear in Table 3.1. Thus, Table 3.1 

is applicable to both basic pole and basic zero configurations. 

v 

Only one more general type of response remained to complete the pro- 
cedures for obtaining the K-jW trace for any given open-loop transfer 
function: namely, the influence upon the trace of poles and zeros con- 

r 

tained in the transfer function in addition to the basic combinations 
already discussed. The influence of these additional poles and zeros 
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BASIC POLE CONFIGURATION 


K- jOJ TRACE 


I. Multiple (real) pole. 

. 1 


1. Even multiplicity: 

a. Trace coincident with 
asymptotes for all 4 10 * 

2. Odd multiplicity: 

a. Trace coincident with 

asymptotes for ico o . 

b. Trace coincident with 

yjKf -axis for : 0 . 


II. Real pole-pole combination 

(with root locus on -axis 

between the poles). 


Coincident with ^V|K| -axis 
from K = 0 to *V K = V 
then rising toward asymptote. 


III. Complex pair of poles. 


Departs -axis at JU) = complex 

part of pole, approaches the 
K- ju> asymptote. 



RULES FOR SHAPE OF K- 4 IP TRACE FOR BASIC POLE CONFIGURATIONS 
OF OPEN-LOOP TRANSFER FUNCTION 
Table 3.1 

waa most apparent graphically by returning to the multiple pole concept. 
Starting with a real pole. of multiplicity three, the influence upon the 
K- Ju> trace of moving one of these three poles along the <T -axis 
while keeping the other two poles fixed was to cause the trace to "sag" 
below the trace obtained for the triple pole. This is shown in Fig. 3.7. 

It was possible to predict the sag intuitively by considering the graphi- 
cal definition of K. Since K is defined (graphically) as the product of 
the distances to poles divided by the product of the distances to zeros, 
it was immediately obvious that, for. a given value of juJ in the lower 
range where the root locus of Fig. 3.7(b) was distinctly separate from 
the asymptote, the value of K must be larger than the value in Fig. 3.7(a). 
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"SAG" EFFECT OF ADDITIONAL POLE UPON K- 4 U) TRACE 

Fig. 3.7 

Consequently, the sag was anticipated. In addition, as K was increased 



to larger values, the poles of Fig. 3.7(b) would appear closer together 
from the point on the root locus where K was measured until, in the limit 
as K approached infinity, the poles would appear as a triple pole. This, 
then, caused the asymptotic return of the K- trace after the sag ef- 
fect;. 

The effect of zeros upon the K- trace was obtained in a similar 

I 

manner, but in order to show the effect upon a single diagram, the number 
of excess poles, n, of the open-loop transfer function had to remain 
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unchanged in order that the K\ -axis remained the same. Therefore 
pole-zero pairs were added. By again considering the graphical defini- 
tion of K, it was apparent that two effects would be evident depending 
upon whether the pole or the zero of the added pair was closer to a 
test point on the root locus. For the case where the zero was closer to 

I 

the root locus, the distance to the pole exceeded the distance to the zero, 
resulting in an overall value of K which was larger than the value before 
the pair was added. Consequently, the sag effect was again produced. 

Where the pole of the pair was closer to the root-locus than the zero, 
the reverse effect resulted due to a reduction in K, which caused the K-juj 
trace to "rise" above the trace obtained without the pole-zero pair. 

These effects are illustrated in Fig. 3.8. 

Based upon the procedures developed for sketching the K- ^ u> trace, 
the following rules were established to be used as a general procedure 
to be followed in obtaining the K-_ju> trace for any given open-loop 
transfer function: 

1. Sketch the root locus plot for the open-loop system. 

2. Group the poles of the transfer function into the 
basic configurations listed in Table 3.1. Poles and 
zeros left over after this basic grouping is performed 
are designated in pairs as ''additional pole-zero pairs". 

3. ■ Poles left over after the pairing operation described 

above are paired together (if there is root locus on 
the (T -axis between them) and designated as "pole-pole 
combinations". 

4. Any real pole remaining is designated as an "unpaired 
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GH - KU*c) 

(a+O) 2 u+ b) (a+c)) 





SAG EFFECTS DUE TO ADDITIONAL POLE AND POLE-ZERO PAIR, K £ 0 AND 4^>0 



rue K(A+d) 

W (*►*>) (a+c) 




RISE EFFECT DUE TO POLE -ZERO PAIR (PORTION OF K- 4 IP TRACE DUE 
TO CIRCULAR ROOT LOCUS NOT SHOWN) 

SAG AND RISE EFFECTS DUE TO ADDITIONAL POLES AND POLE- ZERO PAIRS 

1 

Fig. 3.8 
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pole”. 1 

5. Sketch the K- 4u> trace for uj ■ 0 by utilizing the rules 
outlined in Table 2.1. 

6. Sketch the K- ju> trace for Go ^ o for all basic pole con- 
figurations by utilizing the rules outlined in Table 3.1. 

7. Correct for the rise and sag effects of all additional 
pole-zero pairs (and an unpaired pole, if present). 

8. To achieve greater accuracy if required: 

a. For a second-order system, apply the rules of 

chapter II to obtain an exact K- trace. 

b. For higher-order systems, measure the distances 
to poles and zeros from a few selected points on 
the root locus, then calculate K to achieve greater 
accuracy within any questionable range of K. This 
procedure is closely analagous to the refinement 

of a root locus sketch with a Spi-rule to achieve 

J 

greater accuracy. 

With practice, the K- jco trace for any open-loop transfer function 
can be rapidly and accurately sketched utilizing the rules set down pre- 
viously. Several examples for certain types of transfer functions are 

1 

Steps 3 and 4 may result in zero-zero combinations and an unpaired 
zero if the open-loop transfer function should contain more zeros than 
poles. By the same token, step 2 could result in multiple real zeros. 
However, these eventualities are rare enough in engineering practice that 
they are not specifically mentioned in these rules. Should these situa- 
tions be encountered, all the procedures developed for poles are equally 
applicable to zeros and may be applied in either case. 
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Illustrated in section 3.4. 

3.4 EXAMPLES . The following examples were chosen to illustrate the 
procedures developed for sketching the K- 4 ^ tracet 



a. The root locus for this system was sketched and is 
shown in Fig. 3.9. 

b. Since there are three poles the K- trace must be 



slope equal to sin 60® or 0 . 866 . 

d. The open loop transfer function is composed of a 
double real pole and an "unpaired" pole so that the Re- 
trace would be expected to originate at the origin and sag 



below the asymptote due to the "unpaired" pole. The re- 
sulting K- trace is shown in Fig. 3.10. 
e. To verify the accuracy of the sketch, the gain required 
at Co equal to one was determined from the root locus to 
be equal to, 2. 55. From the K- trace at U> equal one, 



a. The root locus was sketched as in Fig. 3.11. 

b. The open loop transfer function contains a double pole, 
an "additional pole-zero pair" and an "unpaired" pole. 

Thus for © > the trace will originate at the origin 

and will sag below the asymptote due to both the "unpaired" 
pole and the "additional unpaired pole-zero pair". The 
effect of the pole-zero pair, , will be to increase 



1. GH - 






on the A/k- 4 ^ coordinates. 

c. The asymptote on the K- 4 UJ plane will have a 




2. GH “ K 4) 



- 4 - 45 “ 
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the gain required to obtain a particular value of , 

since the zero is nearer to the test point than the pole. 
Thus the K- trace will have a larger sag than was 
found in the previous example. The K- trace is shown 



in Fig. 3.12. 

c. To check the accuracy of the K- jw sketch, the gain 
was compared at Oo equal to one. From the root locus, K 
was equal to 3.3. The K- .(60 sketch gave ^JTT = 1-49, 

or K ■ 3.31. At ti) ■ 3.0 the root locus gain equalled 50 
and the K- \<ji trace gave $J~K ■ 3.67 or K = 49.4. 



3. GH - 



K( A.-m'i 



( 4 . 41 ) ( 4 .* $*(4.+s) 

a. The root locus ia shown in Fig. 3.13. 

b. As in example 2, the open loop transfer function is 
composed of a double pole, an "unpaired" pole, and an 
"additional , pole-zero pair". In this case the "additional 
pole-zero pair" is placed so that the pole is closer to 
the testf point than the zero, which tends to move the K-jco 
trace upward-so far, in fact, that the trace actually 
crosses the asymptote. This upward tendency will be partly 
counteracted by the "unpaired" pple which will cause the 

K- jt*> trace to "sag" near the origin. The resulting K-4 W 
trace is shown in Fig. 3.14. It should be noted that the 
effect of the "additional pole-zero pair" will change as u> 
becomes greater than 4.4. Then, the pole will be farther 
from the test point than the zero, which will tend to make 
the trace sag, but this sag will be small since the test 
point is quite distant from the open-loop poles and zero. 
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ROOT LOCUS PLOT 
FOR , ' 

GH = K- 



(*-+')* (A^3) 
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K- 4 to 
OH 



TRACE FOR 

& 



( 4. v») x (^ 4 8) 










f 



I 

/' 



ROOT LOCUS PLOT 
POR 
GH 



utM 4.+3)(a+5) 







Fig. 3.11 
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4 W A 



> < 



K- 



4W 

GH - 



TRACE FOR 
K(4-M) 

(,A+- (4+3) (A+ S) 




Fig. 3.12 
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ROOT LOCUS PLOT 
POR 

™ . KU+Q 

(A+a) U+*) x (4+S) 
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K- 4OO TRACE FOR 
GH - fcU + O 

U + 2) U* 
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^n(W 



This sag causes the asymptotic return of the trace to the 
upper side of the asymptote. 

3.5 INTERPRETATION OF K- ju> TRACE . Having attained the initial ob- 
jective of determining the K- trace for any open-loop transfer 

function, the results remained to be interpreted with respect to the 
reasons for performing the determination as stated in chapter I. The 
difficulty anticipated in interpretation of the trace was the use of the 
m | k | abscissa for third (and higher) order systems. 

For the purpose of obtaining information about the position of the 
imaginary part of complex roots of the closed-loop system for particular 
values of K, the use of the abscissa was not considered to be ob- 

jectionable. This information is available from the 'ty ( K\ -4^° co- 
ordinates with a simple slide-rule calculation for any desired value of 
K. , 

For the purpose of obtaining sensitivity information over a given 
range of K, however, the abscissa is not so handy. The K- jto 

trace plotted on the ^|\K\ -4^ plane can not be used conveniently to de- 
termine whether the frequency variation is constant, linear, exponential, 
asymptotic, etc., over a given range of K. This handicap may be easily 
overcome, however, due to the fact that the conversion from | R| to 1*1 
can' be readily accomplished with a slide-rule. By choosing a few key 
points on the trace in the > VlK\ plane , the conversion to the K- j(*> 

plane is readily accomplished, and the trace in the latter plane can 
be quickly sketched. Little is sacrificed in speed to obtain the addi- 
tional information provided by the K- 4OO coordinates, and the accuracy 
of the trace is still basically dependent -upon the accuracy of the roof 
locus plot. Speed of plotting must necessarily be inversely proportional 
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to the accuracy achieved, but the majority of the time is devoted to an 
accurate root locus plot and an accurate K- ^ trace in the VTk) 
plane - the conversion of the trace to the K- 4^ plane would involve 
a negligible amount of time in comparison to that amount already spent in 
obtaining the first trace. 

It was concluded, therefore, that the use of the ■yn<r abscissa 
does not present a material disadvantage to the foregoing methods of 
obtaining the K- trace. 
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IV 



CONCLUSION 

The root locus technique is readily extendable into three dimensions. 
Projecting the three-dimensional root locus plot upon the K- plane 
(for second-order systems) or upon the ^TTkT plane (for systems of 

any order) may be quickly and accurately accomplished to obtain the K- 
trace of the three-dimensional locus. Closed-loop system response may 
then be evaluated from the trace with regard to the variation of 

the imaginary part of complex roots with the root locus gain and with re- 
gard to sensitivity over any range of the root locus gain. 

This paper develops and lists the rules for obtaining the K- ^u> 
trace, the utilization of which requires only a knowledge of the pole- 
zero configuration of the open-loop transfer function. Several practi- 
cal examples applying these rules are illustrated in chapters II and III. 

The use of the K- _{u> , trace provides another useful tool for the 
engineer in the design and .analysis of feedback control systems. 
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